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ABSTRACT 

Context, Previous observations with the Infrared Astronomical Satellite and the Infrared Space Observatory, and ongoing observations 
with Spitzer and AKARI have led to the discovery of over 200 debris disks, based on detected mid- and far infrared excess emission, 
indicating warm circumstellar dust. In order to constrain the properties of these systems, e.g., to more accurately determine the dust 
mass, temperature and radial extent, follow-up observations in the submillimetre wavelength region are needed. 
Aims. The /3 Pictoris Moving Group is a nearby stellar association of young (~12Myr) co-moving stars including the classical debris 
disk star /3 Pictoris. Due to their proximity and youth they are excellent targets when searching for submillimetre emission from cold, 
extended, dust components produced by collisions in Kuiper-Belt-like disks. They also allow an age independent study of debris disk 
properties as a function of other stellar parameters. 

Methods. We observed 7 infrared-excess stars in the /? Pictoris Moving Group with the LABOCA bolometer array, operating at a 
central wavelength of 870 pm at the 12-m submillimetre telescope APEX. The main emission at these wavelengths comes from large, 
cold dust grains, which constitute the main part of the total dust mass, and hence, for an optically thin case, make better estimates 
on the total dust mass than earlier infrared observations. Fitting the spectral energy distribution with combined optical and infrared 
photometry gives information on the temperature and radial extent of the disk. 

Results. From our sample,/? Pic, HD 181327, and HD 172555 were detected with at least 3<x certainty, while all others are below 2cr 
and considered non-detections. The image of f3 Pic shows an offset flux density peak located near the south-west extension of the disk, 
similar to the one previously found by SCUBA at the JCMT We present SED fits for detected sources and give an upper limit on the 
dust mass for undetected ones. 

Conclusions. We find a mean fractional dust luminosity f^ ust = 11T0~ 4 at t w 12Myr, which together with recent data at 100 Myr 
suggests an / dust oc r a decline of the emitting dust, with a > 0.8. 

Key words. Stars: circumstellar matter - Stars: individual: /? Pic (HD 39060), HD 172555, HD 181327 - Stars: planetary systems: 
formation - Stars: planetary systems: planetary disks 

1. Introduction tinuously replenishe d by collisions of larger bodies (see e.g. 

iBackman et al.ll995t) . Deeper and more targeted observations by 

The first large and unbiased survey of the infrared (IR) sky, con- the In f rared Space Observatory (ISO), Spitzer, and AKARI, have 

ducted by the Infrared Astronomical Satellite (IRAS), revealed rev i se d and extended t his list of debris disk candidates to en- 

that approximately 15% of nearby main-sequence stars have an compass over 200 stars {M atthews et al] |2(M iRhee et alj|200l: 

excess of 12-100 pm emission, corresponding to a luminosity at |UK ATCD ebris Disk Database 2ods|) 

least 2 x 10" 5 times higher than that expected from a pure stel- Although some of the most nearby debris disks have been 
larphotosphere (|Backman & Paresce| | 1993| ; | Backman & Gillett| imaged in the optical and IR , most are spatiaUy unresolved and 
119871 |Aumannetal-l | 1984i) . Thls indicates that these stars are character i se d solely on fits of the spectral energy distribution 
surrounded by warm circumstellar dust, interpreted as original- (SED) to stellar synthet i c spectra and a handful of IR photom- 
ing in a disk of debris left over from planet formation, which etry data points In order to better constrain the SED, and with 
is being heated by stellar optical and ultraviolet radiation, and this the temperature and radial ext ent of the disk, complemen- 
re-radiating at mid- and far-IR wavelengths. The short orbital tary submillimetr e ( su bmm) photometry is needed. Submm ob- 
hfetime of small dust particles suggests that they are being con- servations probe the Ra yleigh-Jeans tail of the thermal radiation, 
where the measured integrated flux is roughly (sometimes with 



Send offprint requests to: R. Nilsson non-negligible correction) proportional to the temperature and 

'Based on observations with APEX, Llano Chajnantor, Chile (ESO the mass of the dust. Since the most efficient submm emitters are 
programmes 079.F-9307(A) and 079.F-9308(A)). large and cool grains, which also dominate the mass of the disk 
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dZuck erman 2001), observations with in struments like the L arge 
APEX BOlometer CAmera (LABOCA, ISiringo et alll2007h on 
the Atacama Pathfinder Experiment (APEX, Giisten et al. 2006) 
telescope (operating at 870 fim) or SCUBA ^Holland et"alTl 999) 
at the James Clerk Maxwell Telescope (JCMT. |Prestagd ll996) 
(operating at 850 //m) provide good estimators for the total dust 
mass (assuming an optically thin disk at these wavelengths). 

Another advantage of observations in the submm regime is 
the ability to investigate cool dust created in very extended de- 
bris disks or belts located some hundreds of AU from the star 
(akin to the Solar system's Kuiper Belt), compared to the inner 
warm dust at 1-100 AU radius studied in IR surveys. The larger 
disk region probed by submm, compared to IR, makes it feasible 
to potentially resolve nearby disks, and morphologically deter- 
mine oute r disk radii , dynamical interaction with unseen planets, 
etc. (e.g. lWvattll2008l) . even though the angular resolution in the 
submm in general is lower than for shorter wavelengths. 

Young, nearby stars, with a confirmed mid- and far-IR ex- 
cess are prime targets for a submm search f or cold extended 
disks . The /3 Pictoris Moving Group (BPMG, IZuckerman et al.l 
1200 11) is a nearby stellar association of 30 identified young 
(~ 1 2Myr) co-moving member sta rs (IZuckerman & Songl 
2004; Ban-ado v Navascues et al.| [l999). harbouring the perhaps 
most studied of all d ebris disk systems, namely /3 Pictoris 
dSmith & Terrild 11984 . and also e.g., AU Mic. In this paper 
we present a study of 7 main-sequence IR excess stars in the 
BPMG, observed at 870 fan with LABOCA at APEX in Chile. 
The stars were selected from the sample of currently known 
members of BPMG fulfilling the following criteria: (1) located 
at 5 < -50°, thus, southern stars that are not accessible to JCMT, 
but that will be observable with ALMA, and (2) previously de- 
tected IR excess. However, several of the stars previously ob- 
served at 850jum, e.g. AT Mic and AU Mic, had to be excluded 
due to time constraints. By targeting members of the same mov- 
ing group we ensure that any differences in dust properties and 
mass reflect intrinsic, age-indepe ndent, variations, s ince these 
stars are assumed to be coeval (Mentuchet al. 2008). This in- 
vestigation is the precursor of a larger survey of 10-100 Myr 
old stellar associations, which will permit a statistical survey of 
disk properties during crucial time periods of debris disk evolu- 
tion. The ultimate goal is to determine: the incidence of cold dust 
disks around main sequence stars; whether these disks can serve 
as indicators of planetary systems; the physical characteristics, 
e.g. chemical composition and grain size distribution, as a func- 
tion of fundamental stellar parameters such as mass, metallicity 
and age; and the disk lifetimes. 

2. Observations and Data Reduction 

Since May 2007, the 12-m diameter submm telescope APEX, 
operated jointly by Onsala Space Observatory, the Max- 
Planck-Institut fiir Radioastronomie, and the European Southern 
Observatory at an altitude of 5100 m in the Chilean Andes, has 
been equipped with the LABOCA bolometer array. It operates at 
a central wavelength of 870 fim (345 GHz) and a bandwidth of 
150/im (60 GHz), covering a 1 l'A field-of-view with its array of 
295 bolometers; each with a 18'.'6 full- width at half-maximum 
(FWHM) beam. 

Observations were performed between July 27 and August 
5, 2007, employing a spiral pattern scan in order to recover 
fully sampled maps from the under-sampled bolometer array 
(Siri ngo et al.l 12007). Each 7.5 min long scan gives a raw map 
with uniform noise distribution for an area of about 4' around 
the central position of the source. Targets and observation data 



are presented in Table[T] Additional calibration measurements to 
determine the pointing accuracy were made on Jupiter. Uranus 
and Neptune were observed for absolute flux calibration together 
with secondary calibrators B13134, G5.89 and IRAS16293, 
while several skydips in-between target observations yielded the 
correction for atmospheric opacity (with a zenith opacity rang- 
ing between 1.0 and 3.0). 

For data reduction, we used mainly the Java-based software 
MiniCRUSH, a special version o f CRUSH (C omprehensive 
Reduction Utility for SHARC-2, iKovacsl 120081) adapted for 
the APEX bolometers, while IDL and Matlab were used for 
further analysis and calculations. The reduction procedure in 
MiniCRUSH started with a correction for the atmospheric opac- 
ity calculated from observed sky temperature during the skydips, 
with cryostat temperature drifts taken into account. Each tar- 
get scan was corrected with linearly interpolated opacity values 
from its nearest preceding and nearest subsequent skydip scan. 
Thereafter, the target scans were similarly flux calibrated from 
neighbouring calibrator measurements. After these initial cor- 
rections the time-series signal of each bolometer and for each 
scan was reduced by incremental modelling, following the gen- 
eral reduction steps listed below: 

1. 1/f-filtering (over 30 s scales for first iteration, then over 10 
x crossing times) 

2. Corrections for correlated sky noise 

3. Normalisation with mean sky gain (flat-fielding) after flag- 
ging gains outside 0.1-10 

4. Channel weighting 

5. Despiking at 100, 30 and 10cr after flagging channels in the 
array moving too fast or too slow during the spiral pattern 
scan, as well as dead and very noisy channels 

6. Corrections for correlated noise between bolometers sharing 
some parts of the electronics 

7. Normalization to mean bolometer gain (again after flagging 
gains outside 0.1-10) 

8. Noise whitening (in last iteration) 

9. Source mapping by discarding and blanking at 30, 10 and 
3<t (except last iteration), exposure clipping at 10% peak in- 
tegration, and scan weighting via median noise estimation. 

This was iterated 5-6 times to at tain a conver ged solution and 
recover a final map of our target (Kovacs 20081). 

After closer inspection of our calibration-target (point- 
source) scans we noticed non-negligible differences between 
peak fluxes and integrated fluxes hinting at a possibly de-focused 
beam, which should be taken into account when interpreting the 
spatial extent of any detected sources (see sectionHJl. This should 
not have any effect on measured integrated fluxes, however, the 
calibration of data from submm imaging arrays is non-trivial 
since strong, point-source calibrators are not reduced in the same 
way as, e.g. weak extended objects. The differences between the 
reduction procedures have to be quantified i n order to accu rately 
compare source and calibration images (see|Kovacs 20Ql). 

3. Results 

Out of our sample of 7 main-sequence stars in the BPMG, two 
were clearly detected in the reduced 870 /vm maps, while one 
additional star showed what is most likely a source-associated 
submm emission, at a 3>-Acr level. Root-mean-square (RMS) 
noise weighted average maps of the three objects are presented 
in Fig. [T^-c. The colour scale representing the flux level has 
been truncated at zero in order to give cleaner maps. Each im- 
age has been smoothed with a 973-diameter circular Gaussian 
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Target 


R.A. 


Dec. 


Spectral Type 


Distance 


Integration time 




(h m s) 


C ' ") 




(pc) 


(h:min) 


HD 15115 


02 26 16.24 


+06 17 33.2 


F2 


45 


1:45 


HD 39060 (f} Pic) 


05 47 17.09 


-51 03 59.5 


A6V 


19.3 


10:45 


HD 164249 


18 03 03.41 


-51 38 56.4 


F5V 


47 


4:00 


HD 172555 


18 45 26.90 


-64 52 16.5 


A7V 


29 


2:00 


HD 181296 


19 22 51.21 


-54 25 26.2 


AOVn 


48 


1:45 


HD 181327 


19 22 58.94 


-54 32 17.0 


F5/F6V 


51 


2:08 


HD 191089 


20 09 05.21 


-26 13 26.5 


F5V 


54 


1:23 



representing half the FWHM size of the beam, yielding an ef- 
fective image resolution FWHM of 2078. The first solid contour 
corresponds to a 2cr flux level with the following contours at 
increments of la for HD 181327 and HD 172555, and 2cr for 
P Pic, while the dotted contour outlines the l<x level. The po- 
sition of the star as given in SIMBAD astronomical database 
dWenger et alj2000l) corresponds to (Aa, A6) = (0", 0"). The to- 
tal integrated flux density for ft Pic and HD 181327 was found 
by fitting a Gaussian to a defined circular source region, with 
baseline level and gradient correction as measured outside the 
region. For the special case of HD 172555 the integrated flux 
was measured within the rectangles plotted in Fig. [TJ;. Fluxes 
and RMS-noise values <x of the obtained images are presented in 
Table [2] together with the inferred dust temperature, power law 
index of the opacity law, dust mass, and fractional dust luminos- 
ity (see section |4|. For undetected sources an upper (3cr) limit 
on the dust mass is given. HD 15115 was detected at 850 yum 
by [Williams & Andrews! (120061) . who calculated a dust mass of 
3.9MMoon using k^o - 1.7cm 2 g _1 . Considering flux- and tem- 
perature errors, and our use of a slightly higher /(-value, this is 
consistent with our derived upper detection limit, and it probably 
just barely escaped detection in our observations. 

3.1. p Pictohs 

The flux density peak of (3 Pic (Fig. QJ) is 46.7 mJy/beam, cen- 
tered on the position of the star, but with the source just barely 
resolved (with a FWHM of the fitted Gaussian equal to 24". 8), 
causing the dust emission not to appear significantly elongated. 
Previous submm observations with SCUBA at JCMT found a 
margin al extension of the disk with a position angle (PA) of 
32+4° (HollandetaD[T998), which is similar to the PA of the 
circumstellar disk measured from IR a nd optical observations 
(e.g. 30.7+4° from lKalas & Jewitd[l995t) . Also interesting is the 
offset (Aa, AS) = (-26", -22") south-west (SW) 13.1 mJy/beam 
(5<t) flux-peak (labelled I in the figure) with a PA of 36+5° lying 
slightly off the disk plane, at a position roughly consistent with 
the offset (-21", -2 6") "blob" (position labelled A) found by 
iHolland et all (fT998). There is no significant flux at the location 
(labelled B) of the 1200/im peak found bv lLiseau etaf] (120031) . 
The far SW 6<x flux-peak (labelled II) is probably a background 
submm galaxy, but coincidentally has the same PA as our mea- 
sured peak I. Both these features will be discussed further in 
section l4.2.1l 

3.2. HD 181327 

HD 181327 is clearly detected and marginally resolved with a 
fitted FWHM of 31 ".9 (Fig.QJi). The 24.2mJy/beam peak flux 
density is positioned on the star (within the pointing accuracy) 
with a 3<t peak located 22" away at a PA of 107+7° (with the 



uncertainty determined as the angle subtended by the pointing 
accuracy +4" at the peak centre). A Act peak is also found on 
the opposite north-west (NW) extension of the disk, thus, the 
extended submm features are consistent with the orientation of 
the inclined dust ring previously imaged in scattere d light by 
HST/NICMOS (PA of 107+2°; [Schneider et alj|2006l) . In addi- 
tion to the emission along the disk, we find a marginally resolved 
feature orthogonal to the disk, in the southward direction. 

3.3. HD 1 72555 

The intriguing 2-3 cr structure imaged in the vicinity of 
HD 172555 could be a very extended dust distribution associated 
with the star. Considering that a gravitationally bound compan- 
ion, the M0 star CD -64 1208, is located at PA = 56° (almost 
on a 2cr flux-peak, see cross labelled C in Fig. [TJ;) with a pro- 
jected separation of 2000 AU (Zuckerman et al. 2001a), we see 
that the submm flux density distribution displays a striking sym- 
metry relative to the binary system. The integrated flux density 
is 30±10mJy for the SE part of the feature, compared to and 
40±10mJy for the NW side, calculated within the rectangles en- 
closing the brighter features in Fig.QJ. The strongest flux density 
peaks on each side, located at offsets (46",- 14") and (-6", 34"), 
both have a brightness of 19mJy/beam. 

4. Discussion 

4.1. Warm and Cold Debris Disks in the BPMG 

Toget her with the previo usly detected submm emis sion from 
P Pic ( IHolland et al.ll998l) and AU Mic dLiu et al.ll2004l) . our ob- 
servations increase the frequency of detected submm disks in the 
BPMG to almost 17% (5 out of 30 stars). IR observations with 
the Multiband Imaging Photometer for Spitzer (MIPS) of this 
nearby stellar association has suggested a debris disk fraction 
of more than 37%, with 7 members showing a 24 pm excess, at 
least 1 1 with a 70 p m excess, and at least 5 with 160 jum excess 
dRebull et al.ll2008l) . This can be compa red with the -25% disk 
frequency of nearby stars esti mated by iMatthews et aT (2007) 
based on results from IRAS (Backman & Paresce 11991 " and 
Spitzer dRhee et al.l2007l) . Considering the young age (~ 12 Myr) 
of the stars in the BPMG association, the high frequency of 
disks detected both at IR and submm wavelengths is not sur- 
prising. The indication that a large population of cold disks 
may only be detectable at submm wavelenghts (Lestrad e et alj 
2006; IWvatt et al.l l2003) also argues for additional observations 
of members not showing a clear IR excess. The decrease in the 
fraction of disks detected at wavelengths longward of ~ 100 /mi 
is not necessarily an indication of fewer cold disks, but more 
likely a result of lower detection sensitivities, due to decreased 
flux densities at the Rayleigh- Jeans tail of the SED combined 
with lower instrument sensitivities at these wavelengths. Many 
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P Pictoris HD 181327 




R.A. offset, Aa ["] 
(C) 

Fig. 1. Maps of p Pic, HD 181327, and HD 172555 observed at 870 yum with LABOCA at the APEX telescope. The maps have an 
effective resolution of 20'.'8, as shown by the circle. In addition, a scale bar corresponding to 1000 AU at the distance of the star has 
been inserted. The first solid contour represents 2cr flux levels, with the following contours at increments of l<x in (b) and (c), and 
2cr in (a), while the dotted contour outlines the lcr level. The position of the star, as given in SIMBAD, is at (Aa , A6 ) = (0",0"), 
with additional crosses labelled A and B in subfigure (a) representing the "blobs" found by iHolland et al.l d!998l) and iLiseau et al.l 
(2003), respectively. In subfigure (c) the cross labelled C shows the position of the companion star. 



debris disks may have evolved into systems with a very cold 
and low density dust distribution, in which an occasional onset 
of collisio nal dust avalanches from collisions of planetesimals 
or comets dGrigorieva et alj|2007al) . in structures similar to the 
Kuiper-Belt or the inner Oort-cloud, can temporarily increase the 
amount of small particles emitting at submm wavelengths. The 
stochastic nature of such events could explain the wide scatter in 
observed dust masses in disks of the same age and with similar 
spectral types. However, Md us t carries larger uncertanties than 
the flux-derived errors indicated in Table |2l due to uncertanties 
in stellar distance, dust opacity and temperature distribution. A 
more independent value is given by the fractional dust luminos- 
ity, /dust = Liost/L*, which relies only on the flux measurements 
(see section |4~3| >. 

4.2. Spectral Energy Distribution and Disk Properties 

The measured submm flux densities of /? Pic, HD 181327, and 
HD 172555 add an important data point to earlier mid- and far- 



IR photometry, enabling a better determination of the SED of the 
disk, and thereby its temperature and radial extent. The thermal 
radiation emitted from a population of dust particles depends 
on their temperature distribution and opacity. Assuming that the 
opacity index varies with the frequency as a power law, i.e. k v = 
Ko(v/vof, and that the dust grains dominating the flux in the IR 
and submm part of the SED all have the same temperature Td ust , 
then the flux density spectral distribution can be approximated 
with a modified blackbody function of that temperature, as 

2hv 3 ( v V v 3 +^ 

^disk(v) = — - k<)Q = — C. (1) 

( e HST - l) c 2 \ v o/ 

Naturally, the dust grains dominating the submm emission 
would be expected to be colder (located at Kuiper-belt distances) 
than the main IR emitting grains (at asteroid-belt distances), but 
the assumption of equal dust temperature is a generally adopted 
simplification due to limited spatial information. The power law 
exponent, /?, of the opacity law, basically depends on material 
properties and the size distribution of the dust grains. An expo- 
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Table 2. Integrated flux density, root-mean-square noise levels, and derived dust temperature, power law exponent of the opacity 
law, mass, and fractional dust luminosity for the 870-jum observations of main-sequence members of the BPMG. 



Target 


Integrated Flux Density^ F 


RMS Noise, <x 


1 dusU 


B 


Dust Mass^, M dus , 


./dust 




(mjy) 


(mjy/beam) 


(K) 






(10- 4 ) 


/3Pic 


63.6+6.7 


2.4 


89 


0.67 


4.8±0.5 


17 


B Pic (incl. I) 


70.6+8.2 












B Pic (only II) 


15.8±2.4 












HD 181327 


51.7+6.2 


5.0 


70 


0.15 


34±4 


27 


HD 172555 


30+10 (SE^ 


5.2 


320+(10-20) 


0.1+0.1 


10-60 (SE) 


W 




40±10 (NW) 








20-70 (NW) 




HD 15115 




5.1 


62 


0.73£ 


<2.9 


5.81 


HD 164249 




4.0 


70 


< 2.2 


5.9- ; 


HD 181296 




4.8 


150 




< 1.3 


2.4* 


HD 191089 




6.0 


95 




< 3.2 


19~ 



"The integrated flux density was calculated by fitting a Gaussian peak for HD 181327, /3 Pic, and feature II in the /3 Pic map, while a circle of 
radius 47" was used when including feature I. For HD 172555 the flux was measured within the boxes in Fig.QJ, as described i n the text. Estimate d 
errors come from RMS noise calculation in the specific integration region, together with an absolute calibration error of 10% dSiringo e t al. 2009). 

*For detected sources we use the dust tem perature derived from the best fit of the spectral energy distribution. For HD 151 15 the best-fit 
temperature from Williams & Andrews ( 2006) is employed, while temperatures from IR excess measured with IRAS dRhee et al 1 l20Q7h is used 
for the other objects. 

r Upper 3<r limit on the dust mass for undetected sources. 

rf South-east and north-west feature of HD 172555 calculated separately. 

Tor a 15 K disk. 

'From detection at 850 jjm by[wjlliams & Andrews (2006) who found Fgso = 4.9 ± 1.6 mjy and derived a dust mass of 3.8 Mmooh- 
Trom lRebull et alJ J2008h . 
Trom lMoor et ai]( l2006l) . 



nent B ~ would point to grains radiating almost as blackbod- 
ies, while largely unpro cessed interstell ar dust grains have been 
shown to have B « 2 dHildebrandll 19831) . A B of 0.5-1 has been 
found in protoplanetary (T Tauri) disks, indicating amorphous 
or fractal dust grains, or perhaps just gra ins significantly larger 
than those in the interstellar medium dMannings & Emersonl 
119941) . We made ^-minimisation fits to available photometry 
data (weighted with the accuracy of the flux density measure- 
ment) of our three detected objects with the total SED function, 
Fsed(v) = ^ s mr(v)C s tai- + f disk (v)(C d i S k/C), where the stellar pho- 
tosphere F stal is simply approximated by a blackbody function, 
and C s tar and Cdisk are scale factors depending on distance and 
radiating surface area. The best fit for each object is displayed in 
Fig. [2] In Table [731 (available as online material) we have com- 
piled all current mid-IR, far-IR, and submm photometry data of 
the 12 BPMG members observed in the submm so far. At submm 
wavelengths the disk is thought to be optically thin, with the flux 
coming predominantly from larger grains than the ones emitting 
in the IR. An estimate of the dust mass can then be found from 

F D 2 

M dust = v — -, (2) 

KyBy(7dust) 

where F v is the integrated 870 fim flux density, D is the distance 
to the star, and B v can be approximated with the Rayleigh- Jeans 
expression, B v = 2v 2 kT/c 2 , at long wavelengths. We adopt the 
value k — 2 cm 2 g -1 for grain sizes o f ab out 1 mm at A~ 1 mm 
see di scussion in lLiseau et al.1 d2008l) and lMivake & N akagawa 
1993)). Dust temperatures and masses derived from the fits 
are presented in Table |2l where values for undetected sources 
are calculated fr om dust temperatures measured with IRAS 
dRhee et al .1120071) . 

By assuming that the grains are in thermal equilibrium with 
their environment, i.e., that the absorbed stellar radiation energy 



is emitted at the same rate, we can calculate their characteristic 
radial distance from the star as 

The reflectance, A, and emittance, e, of debris disk grains are 
currently uncertain, but assuming amorphous grains without ice 
mantles, we can simplify by approximating them as blackbod- 
ies with a temperature corresponding to a fitted Planck-function. 
This can be justified considering (1) the poro us grains required t o 
explain e.g. polarization maps of AU Mic (Gra ham et alj |2007). 
and outco me from laboratory experiments of dust particle inter- 
actions dWurm & B lum 1998) ; and ( 2) results from e.g. the B Pic 
model by iGrigorieva et alJ d2007bl) showing that any ~ 1 mm 
sized grains at radial distances out to 150 AU would, due to pho- 
tosputtering, be devoid of ice even though located beyond the 
snow-line. 

4.2.1 . y8 Pictoris 

The edge-on debris disk surrounding the nearby A6V star B Pic 
has been intensively studied, with images in optical and near-IR 
light revealing a disk with >1300AU radius, and morphologi- 
cal asymmetries indicating dynamical interaction with planet s 
dKalas & Jewittl 11995b iMouillet et alJ fl997l: iHeap et al.1 boOOl) . 
Further indicati ons of planets, from in ner disk gaps observed 
in the mid-IR (Lagage & Pantinlll994l) was re cently followed 
by ne ar-IR imaging of a possible giant planet (Lagrange et al. 
2009). Sp ectroscopic inv estigations have found silicate emission 
features dChen et al. 2007) and variable redshifted m etallic lines 
suggesting falling evaporating cometary bodies (Beu st et ail 
1994), which could continuously replenish the dust disk and 
expl ain many characteristic s of the near-IR to m illimetre emis- 
sion dLi & Greenberdl998h . lHolland et al.ldl998l) made the first 
submm detection at 850 fim with SCUBA on the JCMT, in which 
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Fig. 2. SEDs of the three detected objects in BPMG. Photometry data-points comes from SIMBAD (B and V mag), 2MASS (J,H and 
K), Spitzer (3.6, 4.5, 5.8, 8.0, 24, 70 and 160 /mi), IRAS (12, 25, 60 and 100 /mi), ISO (60, 90, 120, 150, 170 and 200 fan), JCMT 
(850/mi), SEST (1200 /mi), and our new APEX data (870/mi) marked with a green cross. A stellar blackbody and a disk modified 
blackbody has been fit to the data, yielding T dusl = 89K and p = 0.67 for /3 Pic, T dust = 70K and /3 = 0.15 for HD 181327, and 
r dust = 322 K and p = 0.1 for HD 172555. Subfigure (d) shows the SED of HD 172555 fitted with a simplified stellar (blackbody) 
spectrum, a 320 K modified blackbody dust spectrum, and a second, colder (15 K) dust component. The temperature was derived 
from a simple approximation of blackbody grains at a distance of 1000 AU from the star. The added dust component must have a 
temperature below 20 K in order to match the 870 /mi data without being inconsistent with far-IR data. See section [4~2l for discussion. 



also the SW "blob" was imaged. Our observed submm disk of 
P Pic does no t appe ar as clearly elongated as in the image by 
iHolland et alj dl998l) . but the integrated flux density is similar 
63.6 mJy compared to 58.3 mJy). It is possible that the smaller 
14" FWHM beam of SCUBA at 850/mi allowed the disk to be 
marginally resolved, while the 18" FWHM beam of LABOCA 
did not. Furthermore, the (suspected) slightly out-of-focus beam 
of the present observation would have affected the effective im- 
age resolution. 

The SW flux peak of lHolland etail dl998l) (labelled A in 
Fig. [Tk) is confirmed with a 5cr significance peak, centered at 
the same (within the pointing accuracy of the telescope) posi- 
tion with a PA of 36 ± 5° com pared to 37 + 6°, ju st off the main 
optical disk (PA = 30.7 ±0.7°. lKalas & Jewiti 19951), some 600- 
700 AU from the star. Modelling bv lDent et al] d2000h has shown 
that if this emission component really comes from dust belong- 
ing to the system and is being heated by the star, it must have a 



mass comparable to the rest of the dust disk and would in that 
case also be visible in scattere d optical light, which it is not (see, 
e.g.. iGolimowski et alj|2006h . On the other hand, if the grain s 
are fluffy, as in the comet dust model bv lLi & GreenbergU l998). 
this dust would not have to be visible in scattered light (com- 
pare with the dust ring of e Eri, which is not see n with either 
HST/STIS or near-IR inter ferometry with CHARA: |Proffitt et al.l 
120041: Id i Folco et alj2007l) . In addition to the discussed flux den- 
sity peak, another even stronger emission peak is found at the 
same PA (36 ± 5°) but at a distance of 2000-2500 AU from the 
star. The first interpretation of these bright features would be 
that they are background submm galaxies, but the coinciden- 
tal alignment with features (like warps) in the ob served optical 
and IR disk is puzzling indeed. As discussed in IHolland et al.l 
dl998l) the emission could come from primordial disks around 
younger, low-mass companions to the star. However, in order to 
have avoided detection in the optical and IR such companions 
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must have a mass closer to gas-giant planets. The location of the 
feature observed with SIMBA at 1200yum dLiseau et al.ll2003l) 
does not correspond to any flux exceeding 2<x in our 870/im 
map, however, faint emission can be seen extending southward 
like in the SIMBA and SCUBA images, reaching a 3cr peak on 
the other side of the disk midplane. This latter flux peak is lo- 
cated at a distance comparable to where the most remote dust 
scatte red light has been observed ( ~ 1 400 AU. iLarwoo d & Kalas 
2001) on the SW side of the disk. 

For the source centered flux the best fit (see Fig. |2^) yields 
a dust temperature of 89 K (giving a dust mass of 4.8 MMoon) 
and power law exponent of the opacity law f3 = 0.67. This tem- 
perature is considera bly lower than the 1 30K derived in a sim- 
ilar simple model by iRebull et al.l (120081) but then only consid- 
ering /l<160//m. In addition to our own new submm data, we 
included the 1.2 mm detection bv lLiseau et al. (2003), who also 
arrived at an equally shallow opacity index. The dust in the disk 
of /3 Pic would be expected to have a low /3 as a sign of big 
grains. From Q we get the dust grains characteristic distance 
from the star, which becomes r - 40 AU. This falls in the re- 
gion where dust rin gs have previously been imaged at 17.9/im 
(Wahhai et al. 2003) , and where the Kuiper-belt is located in our 
Solar system (e.g. We issmanlll995l) . However, this is also a re- 
gion which, according to most models, is dust depleted (with 
a peak at ~100-120A U and a relative inner hole inside 50 AU; 
Aug ereau et al.l200ll) . This is an interesting paradox, which goes 
beyond the scope of the present paper, but would have to be ex- 
plained within the frame of a global coherent model of the /3 Pic 
system. 



4.2.2. HD 181327 

The F5/F6V star HD 181327 has, as mentioned earlier, been 
found by HST/NICMOS coronagraphic imaging to have an in- 
clined circumstellar ring of dust. Located at a radial distance 
of 86 AU from the star and being 36 AU wide fechneide r et al.l 
2006) this feature is currently not possible to spatially resolve 
with single-dish telescopes at submm wavelengths since an an- 
gular resolution < 2" would be required. However, Schneider 
et al.'s follow-up observations with HST/ACS revealed a dif- 
fuse scattering halo out to ~500AU, indicating that more ex- 
tended dust surrounds the star. Although our submm image at 
best marginally resolves the source, we do notice an elongation 
in a direction consistent with the PA of the scattering ring, as 
well as southward, almost orthogonal to the disk. In Fig.[2p we 
see that the 870/mi point lies above the best-fit SED, which 
has rd ust = 70K and /3 - 0.15. This suggests that a differ- 
ent and colder population of grains is dominating the excess 
emission in the submm part of the spectrum than that respon- 
sible for the previously dete cted IR excess in IRAS o bservations 
dMannings & Barlowll 19981) and Spitzer IRS spectra dChen et all 
2006). The existence of a n icy Kuiper Belt a round HD 181327 
was recently predicted by IChen et ail d2008l) based on Gemini 
South go-band 18.3-/mi imaging and Spitzer MIPS SED-mode 
observations. Using the higher temperature derived from the fit 
of the modified blackbody function we get a total dust mass of 
34 MMoon, which is remarkably high compared to the few lunar 
masses found in most debris disks. 



4.2.3. HD 172555 

We do not detect any significant submm flux at the position of 
the star, however, the observed features in the region nearby and 



between HD 172555 and its companion star does appear source 
related. The orbital inclination of the A7/M0 binary system is 
unknown, but in order for the observed structure to come from 
dust in bound orbits, it must be locat ed within ~l/3 of the binary 
distance from the primary (see e.g.. lCuntz et al.ll2007l) . thus the 
viewing angle must be fairly high. It is difficult to draw conclu- 
sions on the geometry of the dust distribution, but the similarity 
in the total integrated flux within the NW and SE features, re- 
spectively, and also between the flux densities of their brightest 
features in Fig. QJ could point to dust being produced by colli- 
sions of cometary bodies at the inner rim of a structure compa- 
rable to the proposed dense inner doughnut -shaped part of the 
Oort cloud (lHilislll98U iLevison et alfcOOll) . Although the ob- 
served flux distribution could also originate from e.g., two in- 
teracting background galaxies, the probability of such a chance 
alignment must be considered small. The initial IRAS detection 
of this IR excess so urce revealed an anomalously high dust tem- 
peratu re of 320 K (Oudmaii er et al.|[l992t [Mannings ^; Barlowl 
1 19981) . and Sp itzer IRS spectros copy has shown warm silicate 
dust features dChen et al. 120061) and signs of a recent hyper - 
velocity colli sion between two l arge bodies dLisse et ail 2009). 
Modelling by Wva tt et ail d2007l) suggests that the dust luminos- 
ity is too high to be explained by steady state collisional process- 
ing of planetesimal belts, and that HD 172555 could be under- 
going a transient event. Thus, this system may be experiencing a 
period of heavy bombardment, with frequent collisions of plan- 
etesimal bodies which could explain the high mass of localized 
dust at 1000AU required to explain our observations, and the 
absence of large, cold dust grains at Kuiper-belt distances. The 
main emitting regions are, as mentioned, located at a projected 
distance of ~1000 AU from HD 172555. Consequently, again as- 
suming thermal equilibrium, we can use © to find the expected 
temperature of the dust grains. With the blackbody approxima- 
tion we arrive at 7d us t = 15 K. In addition, comparing the po- 
sition of the 870 fim flux point to the best-fit curve in Fig. [3J:, 
which has 7d us t = 322 K and y3 = 0.10, it can be clearly seen 
that the detected dust must be very cold and that the bulk of 
the submm emission does not come from the same population of 
grains observed in the IR. A second (blackbody) dust component 
with rd ust = 15 K can also be nicely fitted with the SED data (see 
Fig. |2(d)[ >. Although we have basically only one data point sam- 
pling this colder dust, we note that anything above 20 K would 
be difficult to match to the far-IR data, while a significant dust 
population with a temperature below 10 K would be surprising 
at these distances. This range of temperatures has been used to 
estimate the lower and upper limit on the dust mass presented in 
Table [2] More detailed modelling of the extended flux would al- 
low us to determine the geometry and orientation of the emitting 
dust structure. 

For both HD 172555 and HD 181327 additional photometry 
at slightly shorter (e.g. 450/mi with SABOCA at APEX) and 
longer (millimetre with ALMA) wavelengths are important to 
constrain the temperature and nature of the observed colder dust 
component. 



4.3. Fractional Dust Luminosity and Temporal Evolution of 
Debris Disks 

Using our fitted SEDs a nd previous measu r ements by 
Willi ams & Andrews ! d2006l) . IRebull et al.l d2008l) . iMoor et all 
d2006l) . and lLiu et al.l d2004l) . we can make a plot of /dust versus 
spectral type for all currently known debris disks in the BPMG 
(see Fig. |3J. The mean /d ust is 11 • 10~ 4 , with a large standard 
deviation cr = 9 • 10 4 . HD 35850, HD 196982, HD 199143, and 
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HD 358623 all show exceptionally low upper limits, which could 
be due to their special nature; the first being variable, the second 
and third being binary, and the fourth being both. Dynamical in- 
teractions in binary systems would lead to smaller stability zones 
for debris disks and possibly a faster decline in fraction al dust 
luminosity dTrilling et ai1l2007l: iHolman & Wiegertlll999l) . 

A rough comparison of the ~12 Myr BPMG with the upper 
limit, /dust < 2; 10~ 4 , derived for the ~ 100 Myr Pleiades cluster 
dGreaves et al.l l2009) suggests an /dust K Ma usl <* t~ a decline of 
the emitting dust, with a > 0.8. Early results fr om ISO obser- 
vation s of clusters of various ages gave a — 2 { Spangler et al. 
1200 lb . while more recent Spitzer and JC MT data have h inted 
at a more mo dest decline with a = 0.6 dSu et al.l 120061) and 
a = 0.5 - 1 dLiu et al.l 120041) . respect ively. Mode l s of c olli- 
sional dust evolution in debris disks by IWvatt et ail d2007l) are 
consistent with these results, showing an /j ust of 5— 10-10 4 at 
t « 10 Myr, with a falloff as t~ . The relatively large spread of 
/dust also for stars of similar spectral type (and mass) found in 
BPM G (A and F stars) could per haps be explained by the mod- 
els of iKenvon & Bromlevl d2008l) . which show that 10- Myr old 
A-F type stars with a 30-150 AU disk could be just at the steeply 
rising slope of the fractional dust luminosity created by the on- 
set of collisional cascades. The rise from /d ust * 0.1-10 4 to 
100T0~ 4 would come earlier around the more massive stars and 
the stars with more massive disks, perhaps as early as 0.15 Myr 
for a 3 Mq mass star with a 3 MMSN disk, and as late as 10 Myr 
fo r a 1 Mq mass star w i th a 1 /3 MMSN disk (see Figs. 14 and 17 
in lKenvon & Bromlevl (120081) ). 



5. Conclusions 

We summarize the most important findings of our 870 /im ob- 
servations of 7 main-sequence IR-excess stars in the BPMG as 
follows: 

• Out of the 7 stars observed, we made three detections. Two 
of the detected objects have previously never been detected 
at submm wavelengths. Our observations increase the fre- 
quency of detected submm disks in the BPMG to almost 17% 
(5 out of 30 stars). 

• p Pic showed a strong flux density peak centered on the star, 
with two additional weaker flux peaks, both at a PA of ~ 36° 
and radial distances of 600-700 AU and 2000-2500 AU, re- 
spectively. The forme r of these has a position consistent with 
the "blob" imaged bv lHolland eTail dl998l) at 850 yum. Both 
may be background submm galaxies, but surprisingly show 
a close alignment to the observed optical and IR disk plane. 
Simple fitting of a stellar blackbody and a disk modified 
blackbody to the SED data suggests a dust temperature of 
89 K and a low opacity law exponent ji « 0.7. We calculate 
a minimum dust mass using the theoretically derived maxi- 
mum opacity index for 1 mm sized icy grains and arrive at 

4.8M M oon- 

• HD 181327 is clearly detected but at most marginally re- 
solved with an elongation consistent with the previously 
imaged inclined circumstellar dust ring. SED fitting gives 
7/ dust = 70K and /? = 0.15, but with the 870yum flux end- 
ing up above the fit, implying the existence of a population 
of colder dust grains. The calculated minimum dust mass 
-^dust = 34 M Mo on is remarkably high. 

• HD 172555 does not show any significant submm emission 
at the position of the star, however, in the region between the 
star and its lower mass companion a very extended flux den- 
sity distribution is seen. Assuming that the emission comes 



• 


New data 





Previous data 


V 


Upper limits 




- Mean 



Spectral Type 



Fig. 3. Plot of fractional dust luminosities of debris disks 
from fitting SEDs to available flux data of stars in the 
f3 Pic Moving Group, versus stellar spectral type. It shows 
our detections (filled circles): (a) /3 Pic, (b) HD 181327, 

(c) HD 172555, together with previous res ults (open circles): 

(d) HD 15115 dWilliams& Andrews! 120061). (e) HP 16 4249 
dRebull et all 120081). (f) HP 1 8 1 296 dRebull et al J 120081) . (g) 
HD 191089 ( Mooret alJ 120061) . (h) HP 197489 (AU Mic) 
dLiu et al J2004I) . and upper limits (open triangles): (i) HP 35850 
(AF Lep), (j) HP 196982 (AT Mic), (k) HP 199143, and (1) 
HP 358623 (AZ Cap). The upper limits were found from SEP 
fitting of upper (3<x) flux limits (see Table [3] for data and ref- 
erences). For HP 172555 we plotted a fractional dust luminos- 
ity calculated for a second dust component with Tdust = 15 K. 
Mean /j ust (dashed line) and standard deviation (error bar) are 
also shown. 



from dust associated with HP 172555 we derive the thermal 
equilibrium temperature of dust at the projected distances 
(~1000 AU), which becomes 15 K. The integrated flux in the 
observed features lies far above what would be expected for 
a single modified blackbody disk fit to IR photometry data 
in the SEP. Adding a second disk component with a tem- 
perature of 10-20K fits nicely with the observed SEP, and 
agrees with the temperature estimated from dust grain dis- 
tances to the star. The corresponding dust mass would be 
10-60M Mo on and 20-70 M Mo on for the NW and SE feature, 
respectively. 

• We derive the fractional dust luminosity of detected sources, 
which is then plotted versus stellar spectral type. The mean 
fractional dust luminosity, /i ust = 11 -10 -4 measured for de- 
bris disk in the BPMG, agrees with the expected stage of col- 
lisional dust evolution for such ~ 12 Myr systems predicted 
from previous observation s (e.g. ISu et al.l 120061: iLiu et al. 
2004; Spangler et al. 2001) and models ( Ke nvon & Bromley 
2008; IWvatt et al.l 120071) . The large scatter in / dust among 
co-eval stars could be a sign of stochastic collisional dust 
production or a consequence of BMPG's evolutionary stage, 
perhaps just at the onset of collisional du st cascades at 
t « 10 Myr. Comparison with data at 100 Myr dGreaves et al.l 
120091) gives / dust oc r a with a > 0.8. 
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Table .3. Mid-IR, far-IR, and submm photometry of all BPMG members that have been observed at submm wavelengths. A 3cr upper limit on the 
flux density is given for non-detections. 



Star 



Instrument 


A 


Fy 


Referei 




Oum) 


(Jy) 




MIPS 


24 


0.054 + 0.006 


1 


IRAS 


60 


0.473 + 0.052 


1 


ISOPHOT 


60 


0.415 + 0.020 


1 


ISOPHOT 


90 


0.427 + 0.030 


1 


SCUBA 


850 


0.0049 + 0.0016 


2 


LABOCA 


870 


< 0.0153 


3 


MIPS 


24 


0.0790 + 0.00316 


4 


MIPS 


70 


0.0447 + 0.00447 


4 


MIPS 


160 


< 0.077 


4 


IRAS 


12 


0.4900 + 0.0588 


5,6 


IRAS 


25 


0.080 + 0.008 


5,6 


IRAS 


60 


< 0.180 


5,6 


IRAS 


100 


< 1.930 


5,6 


ISOPHOT 


60 


0.049 + 0.014 


5,6 


ISOPHOT 


90 


0.050 + 0.012 


5,6 


SCUBA 


850 


< 0.0045 


7 


MIPS 


24 


7.276 + 0.29104 


4 


MIPS 


70 


12.990 ± 1.299 


4 


MIPS 


160 


3.646 + 0.43752 


4 


IRAS 


12 


2.38 + 0.10 


1 


IRAS 


25 


10.29 + 0.41 


1 


IRAS 


60 


18.92 + 0.76 


1 


IRAS 


100 


9.32 + 0.37 


1 


ISOPHOT 


60 


14.7 ± 0.346 


8 


ISOPHOT 


120 


6.680+ 1.010 


1 


ISOPHOT 


150 


4.390 ± 0.670 


1 


ISOPHOT 


170 


3.807 + 0.143 


1 


ISOPHOT 


200 


2.030 ± 0.340 


1 


SCUBA 


850 


0.0583 ± 0.0065 


9 


LABOCA 


870 


0.0636 ± 0.0067 


3 


SIMBA 


1200 


0.0243 + 0.003 


10 


SEST 


1300 


0.0249 ± 0.0026 


16 


MIPS 


24 


0.076 ± 0.0304 


4 


MIPS 


70 


0.624 ± 0.0624 


4 


MIPS 


160 


0.104 + 0.01248 


4 


IRAS 


12 


0.280 + 0.028 


11 


IRAS 


25 


< 0.158+ 


11 


IRAS 


60 


0.669 ± 0.094 


1 


IRAS 


100 


< 2.63 


11 


ISOPHOT 


60 


0.761 +0.0380 


1 


ISOPHOT 


90 


0.560 ± 0.0390 


1 


LABOCA 


870 


< 0.012 


3 


MIPS 


24 


0.766 + 0.03064 


4 


MIPS 


70 


0.197 + 0.0197 


4 


IRAS 


12 


1.34 + 0.05 


1 


IRAS 


25 


0.91+0.05 


1 


IRAS 


60 


0.241 ± 0.036 


1 


LABOCA 


870 


0.07 ± 0.01 


3 


MIPS 


24 


0.382 + 0.01528 


4 


MIPS 


70 


0.409 ± 0.0409 


4 


MIPS 


160 


0.068 + 0.00816 


4 


IRAS 


12 


0.351+0.018 




IRAS 


25 


0.496 ± 0.025 




IRAS 


60 


0.449 ± 0.040 




ISOPHOT 


60 


0.433 ± 0.022 




ISOPHOT 


90 


0.286 ± 0.021 




LABOCA 


870 


< 0.0144 


3 


MIPS 


24 


0.195 +0.0078 


4 


MIPS 


70 


1.468+0.1468 


4 



HD 15115 



HD 35850 



AF Lep 



HD 39060 



/3Pic 



HD 164249 



HD 172555 



HD 181296 



r\ Tel 



HD 181327 
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Table .3 - Continued 



Star 



Other Name 



Instrument 


A 


1 V 


Refereni 




(urn) 


(Jv) 




MIPS 


160 


0.658 + 0.07896 


4 


IRAS 


12 


164 + 016 

v_/ ■ x v/ r _!_ w« v i u 




IRAS 


25 


286 + 023 




IRAS 


60 


1 93 + 08 




IRAS 


100 


1 69 + 20 




ISOPHOT 


60 


1 730 + 170 




ISOPHOT 


90 


1.410 + 0.140 




ISOPHOT 


170 


736 + 192 




LABOCA 


870 


0517 + 0062 




MIPS 


24 


1856 + 0076 


12 


MIPS 


70 


5443 + 0401 

V/ m^j i r / _i_ v./ . v./ r v./ x 


12 


MIPS 


160 


0.2046 + 0.0446 


12 


IRAS 


25 


387 + 058 




IRAS 


60 


729 + 058 




ISOPHOT 

iu v/i x ivy x 


60 


781 + 081 

\J m 1 \J X _!_ V,/ • V,/ U X 




ISOPHOT 


90 


370 + 038 




LABOCA 


870 


< 0.018 


3 


MIPS 


24 


1160 + 00464 

v./ ■ i i v.' v./ _!_ v./ . v./ v / r\j r 


4 


MIPS 


70 


< 0.018 


4 


SCUBA 


850 


< 0.009 


7 


MIPS 


24 


1430 + 00572 


4 


MIPS 


70 


2050 + 0205 


4 


MIPS 


160 


1680 + 0202 


4 


IRAS 


60 


252 + 043 


1 


SHARC 


350 


0.072 + 0.020 


13, 15 


SCUBA 


450 


085 + 042 


7 


SCUBA 

i-j v_ vj i— / v x 


850 


0144 + 0018 

\j t yJ x i r _!_ \j * \j \j x \j 


7 


MIPS 


24 


0350 + 00014 


4 


MIPS 


70 


< 0.022 


4 


MIPS 


160 


< 0.031 


4 


SCUBA 


850 


< 0.0075 


7 


MIPS 


24 


0.0130 + 0.00052 


4 


MIPS 


70 


< 0.0012 


4 


SCUBA 


850 


< 0.0069 


7 



HD 191089 



HD 196982 AT Mic 



HD 197481 AUMic 



HD 199143 



HD 358623 AZ Cap 
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